Renewable and Sustainable Energy Reviews 32 (2014) 313-327 


Contents lists available at ScienceDirect 


Renewable and Sustainable Energy Reviews 


journal homepage: www.elsevier.com/locate/rser 


Material and manufacturing cost considerations for thermoelectrics Q) cn 


Saniya LeBlanc %P®*1, Shannon K. Yee %1, Matthew L. Scullin €, Chris Dames “, 
Kenneth E. Goodson” 


* The George Washington University, Department of Mechanical & Aerospace Engineering, 801 22nd St. NW, Suite 739, Washington, DC 20052, USA 
> Stanford University, Department of Mechanical Engineering, Building 530, 440 Escondido Mall, Stanford, CA 94305, USA 

€ Georgia Institute of Technology, George W. Woodruff School of Mechanical Engineering, Atlanta, GA 30332, USA 

4 University of California - Berkeley, Department of Mechanical Engineering, 6141 Etcheverry Hall, Berkeley, CA 94720, USA 

€ Alphabet Energy, Inc., 26225 Eden Landing Road, Suite D, Hayward, CA 94545, USA 


ARTICLE INFO ABSTRACT 

Article history: Interest in thermoelectrics for waste-heat recovery and localized cooling has flourished in recent years, but 
Received 6 June 2013 questions about cost and scalability remain unanswered. This work investigates the fabrication costs and 
Received in revised form coupled thermal and electrical transport factors that govern device efficiency and commercial feasibility of the 


18 November 2013 
Accepted 18 December 2013 
Available online 31 January 2014 


most promising thermoelectric materials. For 30 bulk and thin film thermoelectric materials, we quantify the 
tradeoff between efficiency and cost considering electrical and thermal transport at the system level, raw 
material prices, system component costs, and estimated manufacturing costs. This work neglects the cost of 
Keywords: . heat, as appropriate for most waste-heat recovery applications, and applies a power generation cost metric in 
Thermoelectric cost $/W and a cooling operating cost metric in $/kWh. The results indicate material costs are too high for typical 
DE S thermoelectric power generation applications at mean temperatures below 135 °C. Above 275 °C, many bulk 

thermoelectric materials can achieve costs below $1/W. The major barrier to economical thermoelectric power 
generation at these higher temperatures results from system costs for heat exchangers and ceramic plates. For 
cooling applications, we find that several thermoelectric materials can be cost competitive and commercially 


promising. 

© 2014 Elsevier Ltd. All rights reserved. 
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Nomenclature 

A total plate area (m7) 

Aci thermoelectric leg i cross sectional area (m7) 

ATE total thermoelectric material cross sectional area (m?) 

C total capital cost ($) 

Q cost by volume ($/m?) 

c” cost by area ($/m7) 

Cp bulk raw material cost ($/kg) 

Ce price of electricity ($/kWh) 

Cux heat exchanger cost ($/(W/K)) 

Cue manufacturing cost to process bulk material ($/kg) 

Cua manufacturing cost to process material by area ($/m°) 

COP coefficient of performance 

F fill factor 

G generator cost-performance metric ($/W) 

Gmin minimum cost for thermoelectric device ($/W) 

Go cost-performance scaling metric for thermoelectric 
device ($/W) 

H cooler operating cost ($/kWh) 

I electrical current (A) 

k thermal conductivity (W/m-K) 

Ke thermal conductance of cold side heat exchanger (W/ 
K) 

Ky thermal conductance of hot side heat exchanger (W/K) 

Kr total thermal conductance (W/K) 

Kre thermal conductance of thermoelectric legs (W/K) 

Ky thermal conductance of parallel path (W/K) 

E leg length (m) 

Lpr diminishing returns optimum leg length (m) 

Lr ratio of thermal conductivity to heat transfer coeffi- 
cient of heat exchanger (m) 


LCOE levelized cost of electricity ($/kWh) 

m ratio of load resistance to device resistance 

max subscript designating maximum value 

n subscript designating n-type semiconductor 

opt subscript designating optimum value 

p subscript designating p-type semiconductor 

Pen generator electric power output (W) 

Pup electrical power supplied (W) 

Qc heat transferred through cold side heat exchanger (W) 

Qu heat transferred through hot side heat exchanger (W) 

QTE heat transferred through legs (W) 

Qı heat transferred between legs (W) 

r amortization (discount) rate (h7!) 

R electrical resistance (Q) 

R, electrical load resistance (Q) 

S Seebeck coefficient (V/K) 

Spn difference in Seebeck coefficients of p- and n-type legs 
(V/K) 

T temperature (K) 

Tc cold side temperature (K) 

Ty hot side temperature (K) 

Tin mean junction temperature (K) 

Tı hot junction temperature (K) 

T2 cold junction temperature (K) 

AT junction temperature difference (K) 

U heat transfer coefficient of heat exchanger (W/m? K) 

ZT module or material thermoelectric figure-of-merit 

n efficiency 

p density (kg/m?) 

o electrical conductivity (S/m) 


1. Introduction 


Thermoelectric devices are used in power generation and 
cooling applications to either convert heat into electricity or to 
pump heat. To date, thermoelectric technology has been con- 
strained to applications that include vehicle waste-heat recovery 
prototypes, space vehicle power sources, seat coolers, solid-state 
refrigerators, and temperature control in laboratory equipment. 
Thermoelectric conversion has received renewed interest due to 
the development of better-performing materials and their poten- 
tial to improve the efficiency of combustion systems through 
waste-heat recovery [1-3]. Improvements in material performance 
are ongoing [4-6] since conversion efficiencies of typical thermo- 
electric materials remain below 10% [3,7]. 

Thermoelectric devices are attractive because they have no 
moving parts and can be silent, reliable, and versatile. Since 
multiple thermoelectric n- and p-type couples can be connected 
in series, a thermoelectric device can be smaller than a computer 
chip or larger than a solar panel. Nevertheless, considerable 
technical challenges remain. Existing thermoelectric converters 
have lower conversion efficiencies than alternatives such as 
organic Rankine cycles. While some applications are well estab- 
lished, in particular power sources for space vehicles, a variety of 
terrestrial applications - such as high temperature waste-heat 
recovery - are yielding new concerns about reliability and dur- 
ability. Challenges related to sublimation, oxidation, thermal and 
electrical interface degradation, and thermal expansion mismatch 
remain critical for applications of thermoelectric devices [8-15]. 

Thermoelectric technologies face several additional commer- 
cialization challenges. The cost of many thermoelectric materials 


may be prohibitively high, largely due to the use of tellurium and 
germanium in the most common contemporary thermoelectric 
applications [16]. Safety concerns preclude toxic materials such as 
lead [17]. The weight and specific power of thermoelectric devices 
are particularly important for mobile applications such as vehicle 
waste-heat recovery [18], yet few strategies for weight reduction 
have been proposed. While government funding for thermoelec- 
tric technology research and development has expanded signifi- 
cantly in the past decade, the relative lack of private sector 
familiarity with the technology made early stage financing for 
companies comparatively slow to follow [19]. There has been a 
rapid development of materials promising higher efficiencies, in 
many cases through the use of nanostructuring and novel crystal 
structures [1,4,6]. Recent prototypes demonstrated the feasibility 
of thermoelectric devices for widespread, terrestrial, scalable 
applications [20-23]. Multiple start-up companies have received 
recognition and funding while also forming partnerships with 
established academic and industrial research teams [15,18,24—26]. 
The continued development and deployment of scalable thermo- 
electric devices depend on the device cost and energy conversion 
performance [1,27]. 

There is an urgent need for a comprehensive assessment of 
the commercial feasibility of the most promising thermoelectric 
materials while also considering relevant manufacturing and 
system costs. Cost-performance analyses of this type have been 
conducted for other energy technologies such as batteries and 
photovoltaics by considering material extraction costs [28,29]. For 
thermoelectric devices, there have been multiple evaluations of 
the coupled thermal and electrical module performance without 
considering cost [30,31]. The importance of cost was considered by 
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normalizing the material figure of merit by the raw material cost, 
and the costs for specialized applications such as a marine waste 
incinerator have also been determined [17,27,32,33]. A method for 
determining cost using the thermal energy cost and the module 
construction cost for specified material thicknesses has been 
presented [32]. A recent investigation of the cost-efficiency trade- 
offs of thermoelectric generators first determined the geometry 
that optimizes the thermal and electrical performance of the 


Table 1 


device and then calculated the material cost associated with the 
optimized dimensions [27]. The work highlighted the importance 
of optimizing modules in conjunction with the heat source and 
sinks in the overall system. An important extension to the 
optimization approach is made by evaluating the coupled thermal, 
electrical, and economic factors simultaneously [34]. This method 
captures the ability of a moderately efficient device with low 
capital cost to perform better in the market than a highly efficient 


Material identification table. Material type is indicated in the first column and by material name color coding: red for chalcogenides and SiGe, blue for silicides, purple for 
clathrates, green for skutterudites, black for oxides, yellow for half-Heuslers, and grey for others. (For interpretation of the references to color in this table caption, the reader 


is referred to the web version of this paper.) 


Material type ID Material name Manufacturing 
# type 
Chalcogenides and 1 BizTe3 Bulk 
SiGe 2  Bio.s2Sbi.4sTe3 Bulk 
3 Bios52Sby.48Te3 Nanobulk 
4 Bigsaleo.46 Nanowire 
5 (Nao.0283Pbo.945Te0.9733) Nanobulk 
(AginTeosss) 
6 Bi-doped PbSep 9gTeo,o2/PbTe Superlattice 
7 AgPbigSbTe20 Bulk 
8 SiGe Bulk 
9 SiggGe20 Nanobulk 
10 SiGe Nanowire 
Silicides 11. Mg2Sio.ssBio1s Nanobulk 
12 Mg2Sio.6SNo.4 Bulk 
13 -Si Nanobulk 
14 Si Nanowire 
15 MnSii75 Bulk 
16 Mn;5Si28 Nanobulk 
Clathrates 17 BagGaısGezsZN2 Bulk 
18 BagGaı6Ge3z0 Bulk 
19 = Ba7zSr,;Aly6Sizo0 Bulk 
Skutterudites 20 CeFe,Sbj2 Bulk 
21 ~Ybo2Ino.2Co4Sb12 Bulk 
22  Cao.18C03.97Nio.035b12.40 Bulk 
Oxides 23 = (ZNg9sAlo.o2)O Bulk 
24  Caz24Bio.3Nao.3C0409 Bulk 
25  InGaZnO Nanowire 
26 Nag7CoOz-5 Bulk 
Half Heuslers 27  Zro25sHfo.25Tio.sNiSNo.994SĐb0.006 Bulk 
28 = Zro.sHfo.sNio.sPdo.2SNo99Sbo.01 Bulk 
29 Tig. gHfp.2NiSn Bulk 
Other 30  PEDOT:PSS Polymer 


*See Table 2 below for scenario values. 


Material cost Evaluated for low or high temp. scenario Ref. 
($/kg) SSS SS 
ZT in Fopt Lopt Scenario 
material (mm) Temp. 
110 0.74 018 4.53 Low [35,36] 
125 1.05 0.21 441 Low [37 
125 1.52 0.29 3.47 Low [37 
84 0.02 0.07 1.10 Low [38 
81 1.45 0.34 3.01 High [39 
55 1.96 0.02 0.31 Low [40 
84 1.31 0.26 3.59 High [41] 
679 0.30 0.07 2.66 High [42 
371 0.53 0.13 3.39 High [42 
679 0.22 0.06 1.59 Low [43 
6.67 0.67 0.74 29.5 High [44] 
4.04 1.05 ~1 45.5 High [45 
3.09 0.21 0.71 94.5 High [46 
3.09 0.72 0.09 3.38 Low [47 
1.46 0.05 ~1 37.5 Low [48 
1.51 0.07 ~1 35.7 Low [48 
615 0.48 0.13 1.50 High [49 
644 0.36 0.11 1.65 High [50 
1.64 0.09 ~1 38.6 High [51 
37 0.77 0.28 8.34 High [52 
24 0.93 0.31 10.6 High [53 
13 0.77 0.39 17.6 High [54 
2.30 0.08 0.48 58.6 High [55 
30 0.13 0.43 717 High [56 
511 0.07 0.04 1.59 Low [39 
36 0.52 022: 127 High [57 
9.71 1.38 0.49 19.6 High [58 
8.51 0.69 046 15.4 High [59 
10.70 0.41 0.41 16.5 High [60 
0.34 0.01 ~8 101 Low [61 


1: Materials with Fope~1, have high thermal conductivities which suggests that heat exchangers with larger U-values should be used. 
8: The cost of polymers and their manufacturing technique are so different from other thermoelectric materials that a different geometry than conventional TE should be 


considered. 


Table 2 


Scenario parameters, device temperatures, U-values, and heat exchanger (HX) costs. Representative applications are provided for each temperature scenario. 


Applications Thermoelectric Low temperature Solar thermal Automotive exhaust heat Industrial furnace heat 
cooling recovery generator recovery recovery 
Scenario temperature (Cooling) Scenario #1 Scenario #2 Scenario #3 Scenario #4 
Low Medium-Low Medium-High High 
Hot side temperature, 15 100 250 500 800 
Ty (°C) 
Cold side 5 20 20 50 50 
temperature, 
Te (°C) 
Average temperature, 10 60 135 275 425 
Tm (°C) 
Hot side U-value, oo 102 102 120 120 
Up (W/m? K) 
Cold side U-value, oo 105 105 105 105 
Uc (W/m? K) 
HX & Plate Costs - 18.48 18.48 18.48 18.48 


($/(W/K)) 
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device with higher capital cost. While this prior work made 
progress on coupling efficiency and cost, the actual values and 
thus large impacts of manufacturing costs and the costs of key 
system components have not been considered. This is particularly 
important for the newer bulk and thin film materials for which 
manufacturing costs can be high, and new system designs might 
be required to optimize heat exchange. 

The present work provides a detailed investigation of the 
materials, manufacturing, and system costs for 30 leading thermo- 
electric materials [35-61] considering the conversion efficiency. 
We combine these using cost metrics that account for the coupled 
cost and system design implications and trade-offs [34]. These 
metrics enable an assessment of which materials are most pro- 
mising for heat recovery and cooling and incorporate an optimiza- 
tion of material thickness and fill factor. The present work sets 
targets for thermoelectric cooling and power generation. The cost 
of power generation, in $/W, is reported by combining material 
properties, device physics, and raw material costs, all of which are 
reported for seven material classes presented in Table 1 and in 
Appendix F. The analysis is performed for five operating tempera- 
tures shown in Table 2 to reflect the myriad of potential thermo- 
electric applications. The difference in manufacturing costs for 
bulk and thin film materials influences the decision on whether to 
use novel materials. A comparison to existing, competitive tech- 
nologies highlights the cost of current state-of-the-art thermo- 
electric technologies and device performance. 


2. Methods 


This work determines the appropriate device geometry based 
on cost and performance. The fill factor, F, and the leg length, L, of 


a 
F=A,,-/A 

Thermoelectric area, Are = Aen + Aep 
A een ee _ i 
S i p 1 
S ; n-type ! ; P-type ! 
Y I leg 1 I leg I 
& aS Ez== 
a Plate area, A 
ce 
= 
s T; KJ 
S 
S L 
Ba T: n-type | p-type 
3 
n 


R, 
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a thermoelectric device are the dominant design parameters. The 
fill factor is the ratio of the area covered by the active thermo- 
electric material to the plate area, A, as illustrated in Fig. 1. While it 
goes against traditional thermoelectric module architecture, it is 
possible to envision devices with F>1 (see Section 3.1). The 
inverse of the fill factor is known as the thermal concentration 
factor [62]. By reducing the fill factor, less active material is used, 
and the overall cost of the device can be reduced. In many 
applications, the fill factor is fixed to be above some minimum 
by other system constraints (e.g., mechanical robustness, voltage/ 
areal requirements). By reducing the leg length at this fixed fill 
factor, less material is used, and the cost decreases. However, for a 
fixed heat exchanger and fixed F, as the thermoelectric leg length 
is reduced below some optimum, the temperature drop across the 
device decreases, and the output power decreases. This competi- 
tion between power output and material volume suggests the 
existence of an optimum geometry which minimizes cost and 
maximizes power output. Therefore, the most appropriate metric 
for thermoelectric generators should be expressed on a $/W basis 
as 


C 
Pegen’ 


G= a) 
where Pgen is the electric power output of the generator, and C is 
the capital cost of the device. The capital cost includes the 
material, manufacturing and heat exchanger costs which are 
described in detail in Section 2.2. Cost of the heat available for 
conversion is excluded. The contribution of heat exchanger and 
insulating ceramic plate costs are considered in Section 3.2. Both C 
and Pgen are strongly dependent on L and F. A cost-performance 
optimization should be performed on G rather than on C or Peen 
individually [34]. 


P-R/2-IS Qu 


Fig. 1. Diagram of thermoelectric module and the equivalent thermal circuit. The n- and p-type thermoelectric materials are connected electrically in series with a metal 
(shown in orange). The thermoelectric legs are thermally in parallel between electrically insulating ceramic plates (shown in light grey). The analysis is conducted for 
constant hot- and cold-side temperatures T4 and Tc. There are two thermal conductance terms between the junction temperatures T; and T, representing heat transfer 
through the thermoelectric material legs and in the gap between the legs. The internal heat generation term, /7R/2, is the Joule heating component; the Peltier term is ISpnT. 
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A similar argument can be made for thermoelectric coolers. 
In order to reduce the capital cost of the module, shorter leg 
lengths are used. However, this reduces the efficiency of the 
thermoelectric cooler. Since thermoelectric coolers consume elec- 
tricity and provide a cooling service, the cost of electricity is 
important, so a thermoelectric module with a large coefficient of 
performance, COP, is desirable. The COP is the ratio of the cooling 
power removed from the cold side to the electric power con- 
sumed. Electricity is purchased on a $/kW he basis, so the most 
appropriate metric is the operating cost expressed in units of 
$/kW hy, where the subscripts represent electric and thermal 
(cooling) energy, respectively. The capital cost of the thermo- 
electric cooler is amortized over the lifetime and added to the 
operating cost. Operating continuously and purchasing electricity 
hourly over the lifetime allows for a simplification in the amorti- 
zation where the capital cost is multiplied by the amortization 
(discount) rate, r, on an hourly basis (see Appendix B). This allows 
for both capital cost and operating costs to be considered for 
cooling applications. The appropriate operating cost metric is 


Ce C 


H = OPH COP: Ppp 2) 


where Ce represents the price of electricity, and Psy, is the electric 
power supplied to the device to maintain a desired temperature 
difference across the module. The price of electricity varies. In the 
United States, the average industrial user purchases electricity at a 
price of $0.068/kW he, and residential users purchase electricity at 
a price of $0.12/kW he. The average price of electricity is taken as 
$0.098/kW h, in this analysis and is bounded by the industrial and 
residential prices (i.e., the error bars). These values are the average 
U.S. electricity prices reported by the Energy Information Agency 
at the time of article submission. The product of the COP and Psup is 
the thermal cooling power, in kW;y, and yields the thermal energy 
removed in kW hm once divided by r. 


2.1. Device physics 


Both of these cost metrics rely on first modeling the physics of 
the thermoelectric device. Fig. 1 depicts a dual-leg thermoelectric 
module with the legs connected electrically in series and ther- 
mally in parallel. The electrical and thermal transport for this 
geometry has routinely been solved [30,63-66]. The total thermal 
conductance of the thermoelectric module, Ky, is the sum of the 
thermal conductance through the active material, Ky, and a 
parallel thermal leakage conductance, Kı. The heat transferred 
between the legs, Q, is described by Kı which is the sum of 
conduction, convection, and radiation contributions. For most 
applications, conduction is the dominant heat transfer mechanism, 
but radiation and convection can be significant. As a result of these 
three contributions, Kı and Kr are non-linear functions of leg 
length L. A description of how K, is handled in this analysis can 
be found in Appendix D. The heat flowing through the thermo- 
electric legs, Qrg, is augmented by the Joule heating distributed 
throughout the thermoelectric legs and Peltier cooling/heating 
that occur locally at the hot and cold junctions. 

In addition to the thermal conductance of the module, the 
thermal conductance of the hot- and cold-side heat exchangers, Ky 
and Kc, respectively, strongly influence the overall device perfor- 
mance. These conductances are further reduced by any series 
temperature drops due to the electrically insulating ceramic plates 
and module thermal contact resistances. The effect on device 
performance from the heat exchangers has been considered [67]. 
In this analysis, the contribution from the metallization, ceramic 
plates, and heat exchangers is also considered, and details can be 
found in Appendix D. The values of the thermal conductances are 
dependent on the application and heat exchangers (i.e., reservoir 


temperatures and working fluids). The U-value of the heat 
exchangers and the thermal conductance are related to each other 
by the heat transfer area, K=UA, and the values [68] used in this 
analysis along with application reservoir temperatures are sum- 
marized in Table 2. For the scenarios selected in this analysis, the 
appropriate hot side and cold side exchangers have approximately 
the same U-value, so matched heat exchanger simplifications (i.e., 
Ky=Kc=K) are appropriate. However, for exactness these simpli- 
fications are not used in this analysis. Ultimately, these series 
conductances relate the hot- and cold-side reservoir temperatures, 
Ty and Tc, to the hot-side and cold-side junction temperatures, 
Tı and Tp. 


2.1.1. Power generation 

The temperature difference governs the energy that can be 
converted into a thermoelectric generator. For the electrical load 
matched condition, energy balances of the device at the hot and 
cold junctions yield a system of coupled non-linear equations: 


Seal -T S(T —T2)? 


Ky(Ty -T1)—Kr(1 —T2) JR - BR 


0 (3a) 


Soa (Tı ing T2) 
4R 


which can be solved numerically for the exact junction tempera- 
tures. In this expression, Spn=Sp— Sn is the difference in thermo- 
powers of the two legs, and R is the internal electrical resistance of 
the module. Under the approximations that the electrical contact 
resistance, metallization, and shunt resistances are negligible, and 
the thermoelectric legs have the same electrical conductivity 
Op=On=6, the internal electrical resistance is R=4L/(cAF). The 
analysis is conducted assuming the properties are constant within 
each of the five temperature ranges considered, and the properties 
are evaluated at the mean of the reservoir temperatures (i.e., 
Tm=(Ty+Tc)/2). Thermal and electrical contact resistances are not 
considered as these vary widely and are unique to each device; 
they are most important for short leg lengths, typically <1mm. 
These effects can significantly reduce the performance of a device, 
but neglecting these terms provides a lower (optimistic) bound on 
the $/W cost determined for each material. 

The amount of power converted by a thermoelectric generator 


Ky(TyH -T1)+Kc(T2—Tc) = 9, (3b) 


Spn” (T1 -T2 m 


Peen = R (m+ De (4) 
with 

Rı 
m= (5) 


as the ratio of the load resistance, R;, to the internal electrical 
resistance, R, of the module. The thermoelectric generator can be 
operated at maximum power when m=1 (i.e., load matched) or at 
maximum efficiency when m=J/1+ZT». It has been recently 
reported that, due to thermal feedback from the load resistance 
(e.g., increasing the load resistance increases the junction tem- 
peratures difference slightly), m should be slightly greater than 
these nominal values [31,69]. While this is strictly true, the effect 
is a higher order correction, and traditional analysis [30,31] 
suggests that one of these m conditions should be selected. The 
value of m that minimizes the $/W is m=1 because this condition 
maximizes the denominator of Eq. (1). Therefore, the load resis- 
tance is set to exactly equal the internal resistance [34]. 
In this work, ZT is the module figure-of-merit 


= Spn?T 
ai K7R 


(6) 


318 S. LeBlanc et al. / Renewable and Sustainable Energy Reviews 32 (2014) 313-327 


and Tm=(T1+T2)/2 is the average of the junction temperatures. 
Eq. (6) is similar to the material figure of merit ZT=S?oT/k, with 
the material’s thermopower, S, electrical conductivity, o, and 
thermal conductivity, k. The module and material ZT are related 
as discussed in Appendix A. For this manuscript, the material ZT 
will be used to identify the materials in the figures while the 
module ZT is used in the analysis under the simplification that the 
n- and p- type legs have the same thermopower magnitude, same 
electrical conductivity, and same thermal conductivity. 

Using the electric power and the cost of the device, Eq. (1) can 
be evaluated in units of $/W as a function of L for a chosen F as 
shown in Fig. 2a for a doped bismuth telluride example (Material 
ID #2). By selecting a fill factor, an optimum length, Lope exists 
which minimizes G. The heat flowing into the hot side, Qu, and the 
thermal efficiency, 77=Pgen/Qy, are determined as a function of leg 
length as shown in Fig. 2b. Using Eqs. (3)-(5), the electric power 
generated is determined as a function of leg length as shown in 
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Fig. 2. Example of optimum leg length Lopt and cost G minimization for power 
generation using Material ID #2 in Table 1 (bulk Bio.s2Sb1.48Te3) with Fop;=0.21 and 
parameters in the low temperature scenario of Table 2. The leg length that results 
in the minimum cost (a, Lope=4.4 mm) is not the same length that maximizes the 
efficiency (b, Lmax,,=24 mm) or output power (c, Lmax,p=7.3 mm). Note (a) is on a 
linear scale while (b) and (c) are on a log scale. The black and blue y curves in 
(b) are nearly indistinguishable. (For interpretation of the references to color in this 
figure caption, the reader is referred to the web version of this paper.) 


Fig. 2c. The leg length which minimizes G is different from the leg 
length that maximizes efficiency or the leg length that maximizes 
the electric power output. Fig. 2 shows the difference between 
designs optimized for efficiency or power and an alternative 
design that optimizes the ratio of cost to power output. 

The cost, efficiency, and electric power output are functions of 
the load condition m as discussed in Appendix A and as illustrated 
by the different color lines in Fig. 2. In situations where there is 
excess heat (i.e., fuel/heat is free) the maximum power condition 
m=1 is most appropriate and yields the lowest cost. In situations 
where heat is costly, the maximum efficiency condition may be 
more appropriate. Evaluating fuel costs is beyond the scope of this 
work but has been investigated elsewhere [32]. The unmatched 
condition (i.e., m is arbitrary) may also be acceptable depending on 
the application or cost of power electronic conditioning circuits 
necessary to load match. 

Power electronic conditioning circuits match the resistance of 
the load to the resistance of the module in order to maintain either 
the maximum power condition (m=1) or the maximum efficiency 
condition (m=./1+ZTm). The cost-necessity of these circuits is 
evaluated in Fig. 3. The lack of a proper load condition results in an 
opportunity cost for power generation which is plotted as a 
function of ZTm. The opportunity cost is expressed as a percentage 
of the capital cost of the device. The m=1 (horizontal blue line) 
case is the ideal situation where there is no opportunity cost. 
Other arbitrary load conditions (horizontal red lines) represent 
different degrees of opportunity costs. The black curve in Fig. 3 is 
the percent difference between G with m= ./1+ZT, and G with 
m=1; this is the opportunity cost for operating at the maximum 
efficiency condition rather than the load matched condition. The 
black circles represent the position of materials considered in this 
study along the maximum efficiency line. For most materials in 
this study, ZTm is low enough that operating at the maximum 
efficiency or maximum power makes little difference. However, 
at larger ZT,, this difference becomes significant. For a device 
operating at a condition other than m=1, the opportunity cost is 
higher, and the added cost of power electronics should be consid- 
ered. For example, if a material has a ZTm= 1, and throughout normal 
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Fig. 3. Cost associated with non-optimized load conditions. The opportunity costs 
are percentages of the total capital cost of the module. Designing the module to 
operate at the maximum power condition or introducing a power electronics 
conditioning circuit to load match has appreciable cost benefit. The horizontal lines 
represent specific load conditions m=R;/R: (blue, m=1) maximum power condi- 
tion (black, m = /1-+ZTm), maximum efficiency condition (red, m=1.50, 1.75, 2.00) 
unmatched load conditions. The circles represent the practical maximum efficiency 
conditions of the materials investigated in this study. The optimum load is the 
matched load (m=1), which produced the maximum power and the zero 
opportunity cost. (For interpretation of the references to color in this figure 
caption, the reader is referred to the web version of this paper.) 
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operation m fluctuates to 1.5 for a large fraction of generating time, 
then an opportunity cost of 4% of the capital cost is incurred in the 
absences of conditioning circuits. If the conditioning circuit costs less 
than 4% of the capital cost of the device, then the additional power 
delivered warrants the cost of the conditioning circuit. This varies 
with each material/system as they have different capital costs. For 
inexpensive systems the additional system cost of conditioning 
circuits is not justified whereas this cost may be justified for 
expensive systems. 


2.1.2. Levelized cost of energy considerations 

The levelized cost of energy (LCOE) is often used to compare 
the value of different energy generating technologies; a lower 
LCOE value corresponds to a lower cost to generate electricity | 70]. 
This metric is a function of several variables including the over- 
night capital cost, the ongoing fuel costs, and the recurring 
maintenance costs. These costs are spread over the lifetime of 
the system at a specified capacity factor (uptime) and financial 
discount rate. The units of LCOE are typically $/kWh, and 
determining the LCOE value is application-specific because of the 
number (and specificity) of input parameters. The analysis in this 
manuscript only calculates the overnight capital cost of the 
generator (related to G), which is only one input parameter for 
the LCOE value. A simplified LCOE method is presented in 
Appendix C and provides a way of converting the $/W value into 
a $/kW h value using more specific parameters which are beyond 
the scope of this manuscript. 


2.1.3. Thermoelectric cooling 

The device physics of thermoelectric coolers are similar to 
those of generators except that electric power is supplied to 
support a temperature difference rather than being generated. 
It is usually desirable to operate thermoelectric coolers at the 
maximum COP since electricity is continually purchased. For 
thermoelectric cooling the goal is to maintain a specified tem- 
perature difference at the junctions of the module. The hot- and 
cold-side heat exchangers act as additional thermal resistances. 
The cold-side junction Tə needs to be colder than Tc, and the hot 
side junction Tı needs to be hotter than Ty. This over cooling/ 
heating places an extra burden on the thermoelectric module 
design as the cold side junction needs to be cooled below the 
target reservoir temperature. Removing this burden for simplifica- 
tion, ideal hot side and cold side heat exchangers are considered 
here (i.e., V=Uy=Uc=oo). However, the series thermal resistances 
of the ceramic plates and metal shunts are still included, so Ky 
and Kc are finite (see Appendix D). This simplification fixes 
the temperature of the exteriors of the ceramic plates to be the 
reservoir temperatures. The fill factor F is approximately one 
because the device should be covered in as much active thermo- 
electric material as possible in order to pump heat. 

Under these simplifications (U=oo, F=1), the mean junction 
temperature is approximated as the mean reservoir temperature, 
Tm = (Ty+Tc)/2, enabling an explicit solution for the junction 
temperatures by performing an energy balance on the cold side 
junction. This yields a quadratic equation: 
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for which the smallest root is the cold side junction temperature, 
and the hot side junction temperature is determined using 
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Tı x 2Tm—T>2. These junction temperatures are used to determine 
the COP and supplied power P;up: 
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Both of these equations are functions of leg length and can be 
evaluated numerically. Finally, using these equations, the operat- 
ing cost of thermoelectric cooling, H in Eq. (2), can be evaluated. 


2.2. Cost considerations 


In addition to the device physics, both material and manufac- 
turing costs influence the overall cost metric. The cost, C in Eqs. 
(1) and (2), is a function of the bulk raw material cost Cg ($/kg) as 
summarized in Table 1. The manufacturing cost associated with 
processing bulk material Cm, ($/kg), the areal manufacturing cost 
Cua ($/m7), and the cost of both heat exchangers Cyx ($/(W/K)) are 
found in Table F2 in Appendix F. The overall cost C is also 
dependent on the heat exchanger heat transfer coefficient U, the 
density of the active material p, the leg length, and area of the 
ceramic plates, A: 
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The collection of terms in the first expression is a volumetric cost 
C” ($/m?), and the second term is an areal cost C” ($/m?), both of 
which include manufacturing costs. These costs are determined 
based on a conceptual estimating technique [71] since these 
devices are largely unprecedented. This technique projects costs 
for future products which do not exist rather than gathering prices 
for components of a device already in production. Rapid iterations 
on cost estimates completed in tandem with thermoelectric 
product development will prove valuable as devices enter the 
market. 

The cost of the extracted raw, bulk material Cg captures the 
fundamental differences in material costs. The material cost is 
determined from the 2011 price of each element as reported by 
the U.S. Geological Survey; it is based on the worldwide produc- 
tion of the element and represents the average price an industrial 
consumer would pay [72]. Thermoelectric materials are typically 
composed of raw materials which have been processed to a purity 
level of 99% or higher [16]. Typical costs for pure materials are 
presented in Appendix F and can be significantly higher than the 
raw material cost. For instance, aluminum is an abundant material 
with a raw material cost of $2.60/kg. However, the cost of 99.999% 
pure aluminum is approximately $300/kg [73]. While purification 
adds to the material cost [16], the cost of purified materials is 
highly uncertain. It can vary significantly by vendor and is strongly 
influenced by rapidly changing market and processing factors 
[14,59,74,75]. This volatility can convolute the cost-performance 
analysis, so the cost to purify materials is omitted here without 
biasing the results since the purification conditions apply to most 
minerals [29]. A comparison between pure material and raw 
material costs was conducted, and the relative cost differences 
between each material and their subsequent ranking did not 
change when pure material costs were considered. 

Manufacturing costs are divided into two categories based on 
the method of material processing. Processes enacted on the entire 
bulk volume of material such as ball milling and spark plasma 
sintering have costs denoted by Cy,. Other processes like dicing 
and metallization depend on the area of material processed; these 
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costs are indicated by Cya. The cost is calculated by dividing the 
equipment capital cost by the equipment lifetime and throughput. 
The exact cost and throughput values are available in Appendix F. 
In most cases, the cost and throughput of equipment were 
obtained from quotes or estimates provided by equipment ven- 
dors. Throughput is highly dependent on the specific application 
and product, so it is not easily generalized to a variety of materials 
and applications. An activity measure [71] has been assigned to 
the manufacturing equipment considered here which assumes the 
equipment will be in operation year-round. A 20-year equipment 
lifetime and 24h per day, 365 days per year operation were 
assumed. This approach considers the total equipment effective- 
ness [76] and provides the lower bound on manufacturing cost. 
An overall equipment effectiveness which incorporates equipment 
shutdown for issues like repair and maintenance could also be 
considered [76]. The effect on the results presented here are 
minimal for realistic overall equipment effectiveness values. When 
the analysis is conducted assuming equipment operation for 75% 
of the year, the resulting increase in the cost metric G is ~5% or 
lower for most materials. The exception is the superlattice mate- 
rial (Material ID #6) for which the increase is ~ 10%. Additional 
costs of tooling, maintenance, repair, assembly (automated or 
manual), overhead, and labor are not included in this analysis. 

The required throughput of thermoelectric material in a man- 
ufacturing process will depend on both the performance of the 
material and the yield of the process. Fewer legs will be required 
for higher thermopower materials, and less material must be 
processed as yield improves. The amount of active material 
consumed may be different than the amount of active material 
in the final device due to material losses during manufacturing 
processes. For instance, material can be lost during dicing pro- 
cesses, deposited on chamber walls during deposition, or removed 
through etching. The impact of any fractional loss of the material 
on the end cost would be the same as any changes made to the 
equipment effectiveness as described above. A description of how 
to adjust the manufacturing cost for this fractional loss is provided 
in the Appendix F. In order to provide the lower bound for cost, the 
results reported neglect the additional cost of material lost due to 
manufacturing processes. 

Manufacturing equipment costs are approximate and estimated 
for new equipment. Equipment is typically tailored specifically 
for the end-use application, so communication with vendors is 
required for accurate quotes and precise costs. Particularly in early 
stage research, development, and production, tools built in-house 
or used/refurbished equipment may be used to reduce costs. In 
some cases, it may be more cost effective to use a service provider 
or user facility for certain manufacturing steps rather than 
purchasing equipment and performing the process in-house 
[36,77,78]. Equipment and manufacturing costs used herein serve 
as points of reference; real manufacturing costs vary by material, 
process, and application. 

Table 1 summarizes the state-of-the-art materials with pub- 
lished and peer-reviewed thermoelectric properties used in the 
analysis. They are from the following classifications: chalcogen- 
ides, silicides, clathrates, half Heusler alloys, skutterudites, and 
oxides. A novel polymer material, PEDOT:PSS, has recently been 
proposed for thermoelectrics and is also included here [61]. The 
areal manufacturing costs associated with the polymer material 
are unique since a polymer device would likely be made using a 
process like screen printing. A solution printing process can also be 
used for the metallization layer in these cases, eliminating the 
need for dicing. Hence, Cy, for polymer or solution-processed 
materials could be significantly lower than other materials. 

Both bulk and film materials are considered with the latter 
referring to thin film, superlattice, and nanowire materials. The 
raw material cost for each thermoelectric material is also reported 


in Table 1 along with the references for material properties and 
characterization results. Materials are also designated by their 
manufacturing class: bulk, nanobulk, nanowire, superlattice, and 
polymer. Nanobulk materials are ones which were identified by 
the original authors as having nanoscale grain structures. 

The limitations of manufacturing techniques in turn limit the 
achievable range of leg lengths. If the optimum leg length based 
on the cost-performance optimization lies outside the manufac- 
turing limits, entirely new manufacturing methods would have to 
be developed at a significant additional manufacturing cost. Bulk 
materials should be limited to a minimum leg length of 500 um 
as defined by dicing saw limits [79]; however, after executing 
the analysis none of the bulk materials selected for this study 
approached this dicing limit. Furthermore, no leg length limitation 
was placed on nanowire materials since nanowires shorter than 
500 um can be fabricated, and there is no inherent manufacturing 
limitation to making long thermoelectric legs with nanowires. The 
intrinsic, single-nanowire material properties reported in the 
references considered here may be different than a geometry with 
multiple nanowires, but those considerations are case-specific and 
beyond the scope of this analysis. 

Potential thermoelectric power generation applications span a 
large range of operating temperatures. Therefore, the analysis 
was conducted for five operating temperatures representing 
various proposed applications as shown in Table 2. Temperature- 
dependent material properties were extracted from the references 
listed in Table 1, and the material properties at the mean reservoir 
temperature were used in the analysis. 


2.3. Methods summary 


Using the device physics models, material properties, material 
costs, manufacturing costs, and heat exchanger costs, the metrics 
G (Eq. (1)) for power generation and H (Eq. (2)) for thermoelectric 
cooling can be evaluated for a variety of traditional and novel 
thermoelectric materials. In its complete form, G is expressed 
as [34] 
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where the junction temperature difference is evaluated numeri- 
cally by solving a non-linear system of coupled equations (Eqs. 
(3a) and (3b)). A point of diminishing returns in system cost exists 
which designates the optimum fill factor, Fo, [34]. The true 
minimum occurs when F or L is zero. This point is unrealistic as 
no device would exist. Therefore, the diminishing returns point 
nearly minimizes G. The diminishing returns optimum point is 
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By fixing F at this diminishing returns value, Fop, an optimum L 
exists that uniquely minimizes G, designated as Lop. We present 
both the F=F,,; and F=1 scenarios which provide the optimized 
system costs at an optimum Lopt given a fixed F. 

The analysis provides a detailed derivation of a new cost- 
performance metric Go 
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where Lr is the ratio of the material thermal conductivity, k, to 
the heat exchanger thermal conductance, U. Go is a characteristic 
scaling metric for G and accounts for material cost, application 
temperatures, heat exchanger performance, and thermoelectric 
material properties. For a fixed F, when the metal shunt, ceramic 
plate, heat exchanger, and manufacturing costs are excluded, the 
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resulting cost metric Gmin is the minimum material cost for a 
thermoelectric device [34]: 


Gmin = AGoF* (14) 


The fill factor can equal F,,; or some other value as determined by 
device design requirements not considered in this analysis. This 
minimum material cost scales as F°; so exceedingly small fill 
factors will allow for small material costs. For example, for 
Material #2 in the low temperature scenario with a fill factor of 
F=0.01 the minimum material cost is Gmin =$0.0084/W which is 
two orders of magnitude smaller than the value reported in Fig. 5 
with Fopt=0.21. This minimum represents the material cost domi- 
nated regime; however, other regimes where cost is dominated by 
either the areal manufacturing or the heat exchanger also exist 
[34] and have different minimums. 
For cooling, the thermoelectric operating cost H is 
C: CP +CL [~i+ZIm-1 
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where the junction temperature difference and COP are deter- 
mined by solving the quadratic in Eq. (7). In power generation and 
thermoelectric cooling, the thermal conductance of the device, the 
electrical resistance, and the overall cost are complex functions of 
the leg length L and fill factor F. 


3. Results and discussion 


The results specifically address the cost of candidate materials, 
the cost to process those materials into thermoelectric legs for a 
module, and the cost of the heat exchangers making up the full 
thermoelectric device. This enables realistic assessments of cur- 
rent and future materials’ potential for use in power generation 
and cooling applications by linking performance with estimated 
material, manufacturing, and system costs. 


3.1. Material and module manufacturing costs 


The contribution of manufacturing, material, and heat exchan- 
ger costs differ by the manufacturing and material types. Fig. 4 
depicts the cost for processing typical bulk, nanobulk, nanowire, 
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Fig. 4. Example system cost breakdown for various manufacturing classes on a $/W 
basis. These values correspond to the F=Fop; case in Fig. 5a. The low temperature 
case (Scenario #1) was used, and the cost was evaluated for the optimum fill factor 
and leg length (F=F,,;, L=Lopr) for each material. For nanowire materials, the value 
$2500/m? approximating early production microfabrication cost was used in the 
analysis and corresponds to about $1/g for processing silicon. For the superlattice 
material, $1800/m? corresponding to $40/g for processing lead telluride was used. 
Areal manufacturing cost (in green) represents the cost of dicing and the 
metallization of thermoelectric material. The electrically insulating ceramic plate 
and heat exchanger costs are a substantial portion of the overall cost for all 
thermoelectric devices. (For interpretation of the references to color in this figure 
caption, the reader is referred to the web version of this paper.) 


and superlattice thermoelectric materials for power generation. 
The numerical values and cost assumptions are presented in the 
Appendix F. A similar breakdown can be performed for thermo- 
electric cooling, but it is less interesting since the cost of electricity 
dominates. For power generation, the cost of the ceramic plates 
and heat exchangers are a dominant cost for most manufacturing 
classes. Although their contribution to the overall cost is substan- 
tial, heat exchangers can increase the device performance since 
they improve heat transfer to and from the hot and cold reservoirs. 
This analysis suggests that a major reduction in overall system cost 
can be achieved by making inexpensive heat exchangers ( < $1.00/ 
(W/K)) with high U-values ( > 200 W/m? K). Ceramic plates, typi- 
cally used as substrates in planar thermoelectric modules, provide 
the necessary electrical insulation and mechanical support for the 
device and represent an important cost. A significant reduction in 
overall system cost could be achieved by making an inexpensive 
(< $0.10/cm7), thermally conducting (k > 40 W/m K), electrically 
insulating, and structurally suitable substrate material that can 
operate at the given application temperature. 

The advantage of nanobulk materials is evident in Fig. 4 as 
illustrated with a nanostructured bismuth telluride material 
(Material ID #3) compared to the similar bulk material (Material 
ID #2). The nanobulk material has a lower thermal conductivity 
than its bulk counterpart, so it has a lower $/W value. By nanobulk 
processing, a hierarchy of length scales relevant for thermal 
and electrical transport can be engineered to produce a better 
thermoelectric material [80]. These processing techniques are still 
affordable and demonstrate the benefit of a slight increase in 
manufacturing cost that results in an overall reduction in the $/W 
value. This example demonstrates a ~17% overall cost reduction 
on a $/W basis. 

Nanowire materials may require microfabrication processes 
like those used in semiconductor manufacturing. In semiconduc- 
tor device manufacturing, cost typically scales with the number of 
fabrication steps involving lithographic patterning and etching, so 
there is a range in microfabrication costs. The minimum of this 
range corresponds to a high-volume process similar to the current 
state of silicon photovoltaic fabrication. The maximum of this 
range corresponds to multi-mask processes used to make com- 
plementary metal oxide semiconductor (CMOS) devices. In the 
middle cost range are fabrication steps that are used to make 
micro- electro-mechanical systems (MEMS); this is the cost used 
for nanowire manufacturing in this analysis. 

Superlattices are made using molecular beam epitaxy (MBE) 
which requires many hours to deposit a typical thickness of 
material, so the cost associated with achieving sufficiently thick 
thermoelectric legs is high. Since the manufacturing costs asso- 
ciated with thin film materials are significantly higher than those 
for bulk materials, the increase in power output achieved by these 
materials must far outweigh the additional cost. Given the 
manufacturing cost of MBE, little active material is used as 
demonstrated by the small leg length in Table 1. For MBE grown 
materials, bulk material costs are low because little active material 
is used. The cost associated with material deposited on surfaces 
other than the active substrate is not considered here as described 
earlier in Section 2. Appendix F includes details regarding a 
method to adjust the manufacturing cost and compensate for a 
material loss factor. 

Metallization and dicing costs are areal manufacturing costs; 
they are applied to all materials considered herein equally since 
both are required for the traditional approach of assembling 
individual legs into modules. The metallization category denotes 
the deposition of diffusion barriers and metallization to improve 
contact to the electrical shunt. This cost includes both the cost 
of the equipment and an approximate cost of the metal deposited. 
An important exception to traditional metallization and dicing 
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processes is the relatively novel polymer thermoelectric material 
processing. The metallization and formation of individual leg 
structures would likely be done entirely through screen printing 
processes or a similar process commensurate with the unique 
fabrication of polymers. The impact of this exception is evident in 
the final G values computed for this material where the heat 
exchanger cost is neglected, and the system parameters used are 
L=100 mm, F=8, Ty=100 °C. The screen printing approach would 
lead to an estimated polymer thermoelectric module cost of $140/ 
W as opposed to $550/W if traditional metallization and dicing 
costs were applied. 

Large mismatch in costs can yield interesting new results. 
Specifically, the inexpensive material and manufacturing costs 
associated with the polymer and its low thermal conductivity, 
results in an optimum fill factor greater than one. This occurs 
because the cost of the heat exchanger is much larger than the cost 
of the material (i.e., CyxU?>C”k). While a fill factor greater than 
one does not adhere to traditional thermoelectric architecture, it is 
not unimaginable. This is simply the situation where the cross 
sectional area of the legs is larger than the plate area. This is the 
reverse of thermal concentration where heat is actually spread 
laterally through the leg. In this analysis, the same heat exchanger 
was used to compare all materials. This heat exchanger (primarily 
its cost) is not suitable for polymers. This analysis suggests 
polymers require different geometries and heat exchangers than 
traditional thermoelectrics. 


3.2. Thermoelectric power generation 


The cost-performance analysis for thermoelectric power genera- 
tors demonstrates their cost competitiveness when using different 
thermoelectric materials. Table 1 provides an overview of the 
material types, gravimetric costs ($/kg), properties, and optimal 
geometries. In this table, the optimal geometries and ZT are reported 
for a specific temperature application (identified in Table 2 and 
defined by either “low” or “high”) where that material performs best. 
Other temperatures will have a different optimum leg length and fill 
factor for each material. In some cases, particularly nanowire and 
superlattice materials (e.g, Material ID # 6 and 14), only low 
temperature material property data is available, but the optimal ZT 
values may be at higher temperatures. 

There are complex interplays between the thermal transport 
properties, electrical transport properties, and costs that influence 
the thermoelectric leg length. For example, a low thermal con- 
ductivity material allows for shorter legs, which results in lower 
electrical resistances and lower volumetric costs. The consequence 
of high thermal conductivity is starkly illustrated by nanobulk 
silicon (Material ID #13) where the high thermal conductivity 
( >10 W/m K) results in an optimum thermoelectric leg length of 
over 90 mm. A high electrical conductivity permits longer legs, 
which results in larger output power and lower $/W cost. There 
are also complex interplays between manufacturing areal cost, 
heat exchanger cost, heat exchanger thermal conductance, and 
material thermal conductivity that influence the fill factor. For 
example, expensive materials favor small fill factors to reduce 
costs while low thermal conductivity materials favor larger fill 
factors. However, the fill factor is commonly fixed by other system 
constraints, and an optimal thermoelectric leg length is sought. 

This optimum leg length is the length that minimizes the cost 
for a given fill factor (typically F,,;) and is determined for each 
material and temperature scenario. This optimum leg length is not 
the diminishing returns point in Eq. (12) but rather the true leg 
length that minimizes the system cost in Eq. (11) for a fixed fill 
factor. The optimum leg length is used to calculate the values 
reported in Figs. 4 and 5. Three fill factor cases are presented in 
Fig. 5 for each of the temperature scenarios. The first case is for 


F=1, and the second case is evaluated at the optimum fill factor 
F=F py; (Eq. (12)). The last case is a minimum achievable cost 
where F is fixed at the same F,,; value as the second case (i.e., 
Gmin=4GoF opt’, when K, is negligible). This minimum cost is 
reached when the system costs of ceramic plates and heat 
exchangers (Cyx) as well as the metal shunts and manufacturing 
costs (C”) are eliminated, and only the material cost (C”) for the 
thermoelectric device is considered. The gap between this mini- 
mum cost case and the other cases represent the room for 
improvement in manufacturing and system costs. 

While there are many materials that can operate at low 
temperatures for energy scavenging applications, their thermo- 
electric materials costs alone are above $1/W (Fig. 5a). The 
materials with the lowest minimum cost for the low temperature 
scenario, T;,=—60 °C, are predominately bismuth telluride-based 
chalcogenides (Material ID #1-3) although the bismuth telluride 
nanowire material (Material ID #4) has high costs. The presence of 
certain nanowire (Material ID #14) and superlattice (Material ID 
#6) materials at the bottom of the range presented for the 
minimum cost in Fig. 5a highlights the potential of nanostructured 
materials to also be cost competitive. However, the gap between 
the minimum cost points and the other cases indicates that the 
heat exchanger, ceramic plate, and areal manufacturing costs must 
be reduced (to <$1.00/(W/K), <$1.00/cm?, and <$0.01/cm?, 
respectively, see Section 3.4) for most thermoelectric materials 
to be competitive. One option to reduce manufacturing costs of 
nanomaterials may be to use solution-synthesized nanostructures 
since they do not require specialized microfabrication equipment 
and processes. Nanostructured bulk material made with processes 
like ball milling also have lower manufacturing costs than nano- 
wire and superlattice materials. The combined effects of cost and 
performance are responsible for this result; the thermoelectric 
efficiency improvements outweigh the additional costs. 

Certain materials are clearly promising for the mid- to high- 
temperature applications as shown in Fig. 5c. The lead telluride- 
based, half Heusler, skutterudite, and silicide materials designated 
by Material ID #5, 7, 27, 21, and 12, respectively, have consistently 
low costs ( < $6/W). In fact, a factor of two reduction in the system 
costs for these materials would make them highly competitive 
waste-heat recovery solutions as compared to Rankine or organic 
Rankine cycles at $4—5/W, especially when factors such as 
installation costs, reliability, and maintenance are also considered. 
A notably high-cost material type is the oxide category (Material 
ID #23, 24). In spite of the low bulk material cost, the ZT,, is too 
low, largely due to the high thermal conductivity. 

These results should be placed in perspective considering the 
costs of competing electricity generation technologies, especially 
with respect to the heat source temperature [7,81]. The horizontal 
lines in Fig. 5 indicate the cost of competing technologies. When 
the full optimization is performed, Fig. 5a-d shows that current 
thermoelectrics are primarily competitive for mid- to high- 
temperature applications. As shown in Eq. (13), Go is related to 
the inverse of the reservoir temperature difference squared, 
(TH—Tc)~?, which makes the cost of thermoelectric power gen- 
eration sensitive to the reservoir temperature difference. 

Thermoelectric power generating systems are more costly than 
primary power generation sources including coal, natural gas, 
solar, and geothermal. Primary electricity generating sources 
such as coal and natural gas power plants cost less than $3/W 
[82], so it is unlikely that thermoelectric power generation would 
be considered as an alternative to these power sources, especially 
considering the relatively low efficiencies of thermoelectrics 
compared to modern Rankine and Brayton cycles. Through cost 
reductions, renewable energy technologies have recently become 
more competitive. The current cost target for photovoltaics is 
~$1/W, and the photovoltaic module cost is projected to reach 
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$0.50/W [83,84]. Although higher at $3-4/W, the costs of concen- 
trated solar power [84] and geothermal power [82] are still 
competitive and preferable to the costs associated with thermo- 
electric power generation. Key exceptions are cases where com- 
pactness and portability of the generator are critical, and many 
alternatives to thermoelectric power generation are not feasible 
due to resource availability, size, weight, and system complexity. 

Thermoelectrics are especially cost competitive compared to 
other waste-heat recovery technologies. A competitive waste-heat 
recovery technology is the organic Rankine cycle (ORC) system as 
shown in Fig. 5b. The cost for ORC systems is approximately 
$4-5/W [85]. Both ORCs and thermoelectrics can be used to 
recover unused heat. Thermoelectrics can even be complimentary 
to ORCs and scavenge additional heat not recovered by the organic 
Rankine system [86]. The capital costs, system footprint, and 
maintenance costs for organic Rankine systems are considerable, 
so thermoelectric devices may be preferable for waste-heat 
recovery, particularly where reliability is paramount. 

In Fig. 6, the materials can be compared considering two 
parameters of interest: the total areal costs of the module: 


C" module = C"’Lopt +C”, (16) 


and the figure of merit ZTm. The dotted reference lines are drawn 
for comparison to chalcogenide materials: bismuth telluride and 
lead telluride for low and high temperature ranges, respectively. 
These plots provide an image of the current materials landscape 
and can illustrate targets for future materials research. There is a 
paucity of materials in the lower right-hand quadrant which 
indicates new material system modules are needed that are less 
expensive (C’ module < $0.05/cm?) and higher efficiency (ZT > 1) 
than the chalcogenide reference. The clustering of points below 
the chalcogenide comparison material cost point indicates sig- 
nificant research and development work should consider material 


cost while improving thermoelectric efficiency. For example, 
nanowire-based thermoelectric materials requiring microfabrica- 
tion processes are not currently competitive options as shown in 
Fig. 6b and c. Materials which lie in the upper right-hand quadrant 
may be desirable in some instances; military applications such as 
mobile power sources for soldiers or space applications such as 
space vehicle power generators are two examples. 

An important challenge is the selection of materials within a 
given quadrant. When faced with selecting between materials 
which lie in relatively close proximity to one another in Fig. 6, 
there are a variety of application-specific factors that govern the 
selection between, for example, a material having lower cost and 
performance and another having higher cost and performance. 
Such decision-making will likely rely on additional parameters 
related to feasibility such as material stability with respect to time, 
temperature, and compatibility with the operating environment, 
material toxicity, mechanical properties, and compatibility with 
metallization and other module materials, none of which are 
considered in this analysis. 


3.3. Thermoelectric cooling 


Thermoelectric cooling is an application where thermoelectrics 
have performed well in the market. Much of this success can be 
attributed to the room temperature performance of bismuth 
telluride alloys. However, other thermoelectric materials may also 
be attractive alternatives. The cost analysis discussed previously, 
Eq. (2) and (15), was performed on the materials listed in Table 1 
and is presented in Fig. 7. Not all materials are able to achieve the 
targeted AT=T,—T> associated with the example application. 
Eq. (8) shows that if a material does not have a sufficiently large 
ZT for a given AT, then the COP is negative indicating that cooling 
to that AT cannot be achieved. The operating cost in Fig. 7 is 
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expressed for each material at the optimum leg length that 
minimizes Eqs. (2) and (15) (reported in Table F.2). Materials are 
listed by their material identification number and are ordered by 
increasing ZT. There is a clear trend that larger ZTs have lower 
operating cost. This is because the COP is the strongest driving 
parameter, and the purchased electricity dominates the cost. The 
coefficient of performance is directly related to ZT. 

The cost of thermoelectric cooling can be compared against the 
operating cost of absorption refrigeration and air conditioning 
[14,87], and many of the chalcogenides are cost competitive. Material 
ID #6 has the potential to be the least expensive even considering its 
costly MBE process. The primary challenge is making high quality 
superlattices with thickness approaching 130 um, which is the 
optimum length as determined by this analysis. Material ID #3 also 
shows that the improvements in thermoelectric power output and 
associated cost reduction caused by nanostructuring compensate for 
the additional manufacturing expenses. Given the straightforward 
relationship to ZT, an ideal material with ZT ~ oo can have a cost of 
$0.004/kW ha for AT=10 °C. This suggests there is still room for 
material improvements in cooling applications to reduce the cost 
even further. Cost reduction can be achieved by increasing ZT with a 
system capital cost that is recovered by the operating cost savings 
over the lifetime of the device. 


3.4. Thermoelectric device targets summary 


The key targets for thermoelectric devices pertain to system 
parameters and costs. Meeting or exceeding these targets would lead 
to lower cost values than those determined here. As suggested by Go 
in Eq. (13), both large ZT and small volumetric costs C” are desirable, 
but they need not occur simultaneously. Materials with low ZT and 
low volumetric costs can have a lower $/W value than those having 
high ZT and high costs. Hence, there is no unique target for ZT or 
volumetric cost, but approximate targets can be set for current state- 
of-the-art system components. The heat exchanger and ceramic plate 
dominate the overall cost for many devices. For the situation where 
the heat exchanger costs dominate, improving ZT is more important 
than reducing C”. The $1/(W/K) target for heat exchanger cost is 
below the existing minimum cost [68], and U-values above 200 W/ 
m? K are desirable. Ceramic plates with thermal conductivity above 
40 W/m K and cost below $0.10/cm? would reduce thermoelectric 
system costs. Since the areal manufacturing costs contribute mini- 
mally to the overall cost, these costs should remain at or below the 
current value of $0.01/cm? for dicing and metallization of bulk 
materials. Any additional areal processing costs would increase this 
value, and gains in other parameters would need to mitigate the 
additional cost. 


4. Conclusion 


This work applies a new thermoelectrics cost analysis [34] 
which incorporates material properties, device physics, material 
costs, manufacturing costs, and system costs. The analysis and 
resulting cost values provide a tool for thermoelectric device 
research and development. Optimization of thermoelectric leg 
length and fill factor minimizes the ratio of cost to performance, 
as opposed to optimization of efficiency and power output alone. 
This work makes progress on evaluating the status of thermo- 
electric materials relative to competitive power generation tech- 
nologies for a range of application temperatures. Novel nanowire 
and superlattice materials have the potential to have a low $/W 
value if improvements in ZT are made above what is reported, but 
with the currently reported values they are not competitive in the 
near-term due to the large costs associated with microfabrication/ 
MBE manufacturing techniques. When applied to several 


traditional and new thermoelectric material classes, the analysis 
demonstrates the paucity of materials which are both higher 
performing and lower cost than standard chalcogenides. However, 
there are multiple thermoelectric materials which could yield 
power generation systems that are affordable and competitive 
with other forms of power generation; the key challenges for these 
materials will be in both engineering devices around them and 
scaling their production. The utility of the cost metrics presented 
here will evolve as factors such as resource availability and energy 
costs change over time demonstrating the need to revisit this 
analysis periodically. The realization of commercial thermoelectric 
devices depends on the simultaneous improvement and optimiza- 
tion of material properties, system design, system costs, and 
material and manufacturing costs. 
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